Adaptive isogeometric finite element analysis of steady‐state groundwater flow          下载免费PDF全文

Yared W. Bekele  Trond Kvamsdal  Arne M. Kvarving  Steinar Nordal 《国际地质力学数值与分析法杂志》2016,40(5):738-765

Numerical challenges occur in the simulation of groundwater flow problems because of complex boundary conditions, varying material properties, presence of sources or sinks in the flow domain, or a combination of these. In this paper, we apply adaptive isogeometric finite element analysis using locally refined (LR) B‐splines to address these types of problems. The fundamentals behind isogeometric analysis and LR B‐splines are briefly presented. Galerkin's method is applied to the standard weak formulation of the governing equation to derive the linear system of equations. A posteriori error estimates are calculated to identify which B‐splines should be locally refined. The error estimates are calculated based on recovery of the L₂‐projected solution. The adaptive analysis method is first illustrated by performing simulation of benchmark problems with analytical solutions. Numerical applications to two‐dimensional groundwater flow problems are then presented. The problems studied are flow around an impervious corner, flow around a cutoff wall, and flow in a heterogeneous medium. The convergence rates obtained with adaptive analysis using local refinement were, in general, observed to be of optimal order in contrast to simulations with uniform refinement. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

Decadal fingerprints of freshwater discharge around Greenland in a multi-model ensemble   总被引：2，自引：2，他引：0  

Didier Swingedouw  Christian B. Rodehacke  Erik Behrens  Matthew Menary  Steffen M. Olsen  Yongqi Gao  Uwe Mikolajewicz  Juliette Mignot  Arne Biastoch 《Climate Dynamics》2013,41(3-4):695-720

The recent increase in the rate of the Greenland ice sheet melting has raised with urgency the question of the impact of such a melting on the climate. As former model projections, based on a coarse representation of the melting, show very different sensitivity to this melting, it seems necessary to consider a multi-model ensemble to tackle this question. Here we use five coupled climate models and one ocean-only model to evaluate the impact of 0.1 Sv (1 Sv = 10⁶ m³/s) of freshwater equally distributed around the coast of Greenland during the historical era 1965–2004. The ocean-only model helps to discriminate between oceanic and coupled responses. In this idealized framework, we find similar fingerprints in the fourth decade of hosing among the models, with a general weakening of the Atlantic Meridional Overturning Circulation (AMOC). Initially, the additional freshwater spreads along the main currents of the subpolar gyre. Part of the anomaly crosses the Atlantic eastward and enters into the Canary Current constituting a freshwater leakage tapping the subpolar gyre system. As a consequence, we show that the AMOC weakening is smaller if the leakage is larger. We argue that the magnitude of the freshwater leakage is related to the asymmetry between the subpolar-subtropical gyres in the control simulations, which may ultimately be a primary cause for the diversity of AMOC responses to the hosing in the multi-model ensemble. Another important fingerprint concerns a warming in the Nordic Seas in response to the re-emergence of Atlantic subsurface waters capped by the freshwater in the subpolar gyre. This subsurface heat anomaly reaches the Arctic where it emerges and induces a positive upper ocean salinity anomaly by introducing more Atlantic waters. We found similar climatic impacts in all the coupled ocean–atmosphere models with an atmospheric cooling of the North Atlantic except in the region around the Nordic Seas and a slight warming south of the equator in the Atlantic. This meridional gradient of temperature is associated with a southward shift of the tropical rains. The free surface models also show similar sea-level fingerprints notably with a comma-shape of high sea-level rise following the Canary Current.  相似文献   

The genesis of Pb‐Zn sulphide occurrences in the Lower Devonian Buchan Group,Victoria     

D. C. Arne  P. W. Cromie  J. A. Webb  J. R. Richards 《Australian Journal of Earth Sciences》2013,60(2):75-90

Six minor sulphide occurrences hosted by the Lower Devonian Buchan Group have been investigated. Sulphide minerals and associated phases are hosted by both dolostone and limestone lithologies along stylolitized bedding planes, cross‐cutting fractures, low angle minor faults and in cavities. Mineralization was closely associated with minor structures of inferred Tabberabberan age (Middle Devonian), which it appears to have post‐dated, but was largely strata‐bound in nature. The mineralogy of the occurrences is simple and characterized by the following generalized paragenesis which reflects the increasing oxidation state and pH of the mineralizing fluids: pyrite (pseudomorphous after marcasite)‐galena ± sphalerite ± pyrite‐dolomite ± barite‐calcite ± fluorite ± dolomite.

The sulphur isotope composition of sulphide minerals varies from ‐32.1 to +4.1‰, with iron and base metal sulphide minerals forming two distinct populations around ‐25‰ and 0‰, respectively. A single barite sample gives a sulphur isotope composition of +22.4‰, which is similar to that estimated for Early Devonian seawater. Fluid inclusions in fluorite and calcite homogenize at temperatures in the range 160 to 212° C and have average salinities of approximately 10 wt% NaCI eq. Sphalerite contains up to 1.81 wt% iron which correlates with colour, and up to 1.43 wt% cadmium. The Pb isotopic pattern of galena suggests a source region with U/Pb(= μ) lower than the crustal average, and a high Th/U.

A genetic hypothesis is proposed which involves the circulation of saline fluids through the Snowy River Volcanics, which directly underlie the Buchan Group, during or at some time after the Tabberabberan Orogeny. Although the Buchan occurrences show features characteristic of both Mississippi Valley‐type and stratiform ore deposits, they are most directly comparable to the epigenetic zones of Irish carbonate‐hosted base‐metal deposits. However, Pb‐Zn sulphide mineralization at Buchan appears to have been associated with minor compressional structures, suggesting that a simple correlation with the Irish deposits is not directly applicable.  相似文献   

Consistent joint elastic‐electrical differential effective‐medium modelling of compacting reservoir sandstones     

Erling Hugo Jensen  Leiv‐J. Gelius  Tor Arne Johansen  Zhong Wang 《Geophysical Prospecting》2013,61(4):788-802

  相似文献   

Discussion and Reply Timing relationships between felsic magmatism and mineralisation in the central Victorian gold province,southeast Australia     

D. H. Taylor  R. A. Cayley  A. H. M. VandenBerg  V. J. Morand  F. P. Bierlein  D. C. Arne 《Australian Journal of Earth Sciences》2013,60(3):587-591

  相似文献   

Flux comparison of Eulerian and Lagrangian estimates of Agulhas leakage: A case study using a numerical model     

Erik van Sebille  Peter Jan van Leeuwen  Arne Biastoch  Wilhelmus P.M. de Ruijter 《Deep Sea Research Part I: Oceanographic Research Papers》2010,57(3):319-327

Estimating the magnitude of Agulhas leakage, the volume flux of water from the Indian to the Atlantic Ocean, is difficult because of the presence of other circulation systems in the Agulhas region. Indian Ocean water in the Atlantic Ocean is vigorously mixed and diluted in the Cape Basin. Eulerian integration methods, where the velocity field perpendicular to a section is integrated to yield a flux, have to be calibrated so that only the flux by Agulhas leakage is sampled. Two Eulerian methods for estimating the magnitude of Agulhas leakage are tested within a high-resolution two-way nested model with the goal to devise a mooring-based measurement strategy. At the GoodHope line, a section halfway through the Cape Basin, the integrated velocity perpendicular to that line is compared to the magnitude of Agulhas leakage as determined from the transport carried by numerical Lagrangian floats. In the first method, integration is limited to the flux of water warmer and more saline than specific threshold values. These threshold values are determined by maximizing the correlation with the float-determined time series. By using the threshold values, approximately half of the leakage can directly be measured. The total amount of Agulhas leakage can be estimated using a linear regression, within a 90% confidence band of 12 Sv. In the second method, a subregion of the GoodHope line is sought so that integration over that subregion yields an Eulerian flux as close to the float-determined leakage as possible. It appears that when integration is limited within the model to the upper 300 m of the water column within 900 km of the African coast the time series have the smallest root-mean-square difference. This method yields a root-mean-square error of only 5.2 Sv but the 90% confidence band of the estimate is 20 Sv. It is concluded that the optimum thermohaline threshold method leads to more accurate estimates even though the directly measured transport is a factor of two lower than the actual magnitude of Agulhas leakage in this model.  相似文献   

The crustal evolution of South America from a zircon Hf‐isotope perspective          下载免费PDF全文

Carita Augustsson  Arne P. Willner  Tobias Rüsing  Hans Niemeyer  Axel Gerdes  Christopher J. Adams  Hubert Miller 《地学学报》2016,28(2):128-137

Hf‐isotope data of >1100 detrital zircon grains from the Palaeozoic, south‐central Andean Gondwana margin record the complete crustal evolution of South America, which was the predominant source. The oldest grains, with crustal residence ages of 3.8–4.0 Ga, are consistent with complete recycling of existing continental crust around 4 Ga. We confirm three major Archaean, Palaeoproterozoic (Transamazonian) and late Mesoproterozoic to early Neoproterozoic crust‐addition phases as well as six igneous phases during Proterozoic to Palaeozoic time involving mixing of juvenile and crustally reworked material. A late Mesoproterozoic to early Neoproterozoic, Grenville‐age igneous belt can be postulated along the palaeo‐margin of South America. This belt was the basement for later magmatic arcs and accreted allochthonous microcontinents as recorded by similar crustal residence ages. Crustal reworking likely dominated over juvenile addition during the Palaeozoic era, and Proterozoic and Archaean zircon was mainly crustally reworked from the eroding, thickened Ordovician Famatinian arc.  相似文献   

Airborne Gravimetry Survey for the Marine Area of the United Arab Emirates     

Rene Forsbergbi  Arne Vestergaard Olesen  Adel Alshamsi  Arne Gidskehaug  Sahrum Ses  Majid Kadir 《Marine Geodesy》2013,36(3):221-232

The Military Survey Department (MSD) of the United Arab Emirates (UAE) undertook an airborne gravity survey project for the marine area of the country in 2009, especially to strengthen the marine and coastal geoid in the near-shore regions. For the airborne gravity survey, 5 km spacing coast-parallel flight lines were planned and surveyed. These lines were supplemented by cross-lines in order to assess the quality of the airborne gravity surveys. The flight lines were extended 10 km, spacing lines further offshore. A Beech King Air 350 aircraft was used for the surveys, collecting data at a typical flight speed of 170 knots and a typical flight elevation of 900–1500 m, depending on weather conditions and topography. Gravity was measured with a ZLS-modified LaCoste and Romberg gravimeter (S-99), augmented with a Honeywell strap-down inertial navigation system unit. The estimated accuracy for the airborne gravity data is better than 2.0 mGal r.m.s., as judged from the airborne track crossovers. The new airborne gravimetry data changed the UAE coastal geoid by up to 30 cm in some regions, highlighting the importance of airborne gravity coastal surveys.  相似文献   

Modeling the role of rainfall patterns in seasonal malaria transmission     

Arne Bomblies 《Climatic change》2012,112(3-4):673-685

Seasonal total precipitation is well known to affect malaria transmission because Anopheles mosquitoes depend on standing water for breeding habitat. However, the within-season temporal pattern of the rainfall influences persistence of standing water and thus rainfall patterns can also affect mosquito population dynamics in water-limited environments. Here, using a numerical simulation, I show that intraseasonal rainfall pattern accounts for 39% of the variance in simulated mosquito abundance in a Niger Sahel village where malaria is endemic but highly seasonal. I apply a field validated coupled hydrology and entomology model. Using synthetic rainfall time series generated using a stationary first-order Markov Chain model, I hold all variables except hourly rainfall constant, thus isolating the contribution of rainfall pattern to variance in mosquito abundance. I further show the utility of hydrology modeling using topography to assess precipitation effects by analyzing collected water. Time-integrated surface area of pools explains 70% of the variance in simulated mosquito abundance from a mechanistic model, and time-integrated surface area of pools persisting longer than 7?days explains 82% of the variance. Correlations using the hydrology model output explain more variance in mosquito abundance than the 60% from rainfall totals. I extend this analysis to investigate the impacts of this effect on malaria vector mosquito populations under climate shift scenarios, holding all climate variables except precipitation constant. In these scenarios, rainfall mean and variance change with climatic change, and the modeling approach evaluates the impact of non-stationarity in rainfall and the associated rainfall patterns on expected mosquito activity.  相似文献   

Long-term ice sheet–climate interactions under anthropogenic greenhouse forcing simulated with a complex Earth System Model     

Miren Vizcaíno  Uwe Mikolajewicz  Matthias Gröger  Ernst Maier-Reimer  Guy Schurgers  Arne M. E. Winguth 《Climate Dynamics》2008,31(6):665-690

Several multi-century and multi-millennia simulations have been performed with a complex Earth System Model (ESM) for different anthropogenic climate change scenarios in order to study the long-term evolution of sea level and the impact of ice sheet changes on the climate system. The core of the ESM is a coupled coarse-resolution Atmosphere–Ocean General Circulation Model (AOGCM). Ocean biogeochemistry, land vegetation and ice sheets are included as components of the ESM. The Greenland Ice Sheet (GrIS) decays in all simulations, while the Antarctic ice sheet contributes negatively to sea level rise, due to enhanced storage of water caused by larger snowfall rates. Freshwater flux increases from Greenland are one order of magnitude smaller than total freshwater flux increases into the North Atlantic basin (the sum of the contribution from changes in precipitation, evaporation, run-off and Greenland meltwater) and do not play an important role in changes in the strength of the North Atlantic Meridional Overturning Circulation (NAMOC). The regional climate change associated with weakening/collapse of the NAMOC drastically reduces the decay rate of the GrIS. The dynamical changes due to GrIS topography modification driven by mass balance changes act first as a negative feedback for the decay of the ice sheet, but accelerate the decay at a later stage. The increase of surface temperature due to reduced topographic heights causes a strong acceleration of the decay of the ice sheet in the long term. Other feedbacks between ice sheet and atmosphere are not important for the mass balance of the GrIS until it is reduced to 3/4 of the original size. From then, the reduction in the albedo of Greenland strongly accelerates the decay of the ice sheet.  相似文献